Theory of electric dipole spin resonance in quantum dots: 
Mean field theory with Gaussian fluctuations and beyond 
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Very recently, the electric dipole spin resonance (EDSR) of single electrons in quantum dots 
was discovered by three independent experimental groups. Remarkably, these observations revealed 
three different mechanisms of EDSR: coupling of electron spin to its momentum (spin-orbit), to 
the operator of its position (inhomogeneous Zeeman coupling), and to the hyperfine Overhauser 
field of nuclear spins. In this paper, I present a unified microscopic theory of these resonances in 
quantum dots. A mean field theory, derived for all three mechanisms and based on retaining only 
two-spin correlators, justifies applying macroscopic description of nuclear polarization to the EDSR 
theory. In the framework of the mean field theory, a fundamental difference in the time dependence 
of EDSR inherent of these mechanisms is revealed; it changes from the Rabi-type oscillations to a 
nearly monotonic growth. The theory provides a regular procedure to account for the higher nuclear- 
spin correlators that become of importance for a wider time span and can change the asymptotic 
behavior of EDSR. It also allows revealing the effect of electron spin dynamics on the effective 
coupling between nuclear spins. 

PACS numbers: 71.70.Ej, 73. 21. La, 73.63.Kv, 76.20.-|-q 



I. INTRODUCTION 

One of the principal avenues of semiconductor spin- 
tronics is based on the electrical manipulation of electron 
spins in single and double quantum dots that are envi- 
sioned as prospective blocks for quantum computation J^ 
There are two major aspects of this problem. First is 
based in electrical control of spin populations by means 
of the Coulomb and Pauli blockade when electric current 
passing through the dot populates or depopulates specific 
spin states.—'^ Second aspect is based on direct electrical 
operation of electron spin by means of the electric dipole 
spin resonance (EDSR).— i^ Recently, it was achieved in 
quantum dots by three different mechanisms. ^'^'^ 

Nowack et al^ observed Rabi oscillations driven by the 
electric field generated by ac gate voltage and coupled to 
electron spin via spin-orbit (SO) interaction. Laird et 
al^ reported hyperfine-mediated gate-driven EDSR re- 
markable for a nearly monotonic increase of spin polar- 
ization rather than Rabi oscillations of it. The under- 
lying mechanism is spacial inhomogeneity of the Over- 
hauser field acting on electron spin. It is physically al- 
lied to the EDSR mediated by a spatially dependent Zee- 
man Hamiltonian, due to the spatial dependence of either 
the externaliSiiiiii^ or exchange-S'^^ magnetic field or of 
the Lande g-iactoii^ The SO mechanism may dominate 
in strong external magnetic fields B but is suppressed 
in weak fields because the Kramers' theorem requires 
(in confined geometries) the Rabi frequency to vanish 
linearly in B as B -^ q ^15,16^i7,18 -^^j^iigt the hyperfine 
mechanism survives in the B — > limit due to the bro- 
ken time-inversion symmetry and therefore dominates in 
weak magnetic fieldsi^ More recently, Pioro-Ladriere et 
al^ achieved EDSR in a double dot by employing a spa- 



tially inhomogeneous (slanting) stray field of a micro- 
magnet and proved high efficiency of this approach. The 
traditional electron spin resonance (ESR) driven by an ac 
magnetic field has also been achieved in quantum dots^^ 
however EDSR promises higher efficiency and provides 
easier access to individual spins at nanoscale. Golovach 
et al. estimated relative intensities of ESR and SO- 
mediated EDSR in quantum dots and concluded that 
electron spin can be operated at a timescale of 10 nsj^i 

The research described in this paper was inspired by 
the observation of single-electron EDSR in GaAs quan- 
tum dots.^'*'^ Electron spin dynamics in such dots is a 
challenging problem because of the hyperfine coupling of 
the electron spin to nuclear spin bath. As applied to 
spin relaxation in quantum dots, a mean-field theory ap- 
proach advanced by Merkulov et al^, Khaetskii et al.^^ 
and Erlingsson and Nazaro\^ and based on a large num- 
ber of nuclear spins in the dots, typically N w 10^ — 10^, 
proved rather successful. However, some important as- 
pects of electron spin relaxation and dephasing cannot be 
described in the framework of the mean-field approach, 
and dynamics of a single electron spin in a nuclear spin 
environment developed recently into a rather extensive 
field, e.g., see papers ^illgilg^iilail and references 
therein. 

Applying the mean field approach to EDSR is highly 
attractive because it simplifies the problem tremen- 
dously. This approach is based on (i) statistical ar- 
guments (large number of nuclear spins) and (ii) slow 
nuclear-spin dynamics. However, large electron spin- flip 
frequency that keeps nuclear dynamics slow nearly van- 
ishes in the rotating fram o^°'^^ under the EDSR condi- 
tions. Therefore, the criteria of the applicability of the 
mean field approach are far from obvious. In particular, 



the possibility of accelerating spin relaxation in the nu- 
clear bath in the EDSR regime must be examined. We 
come back to this problem in Sec.lVIIwhere it is discussed 
from the standpoint of the results derived in the paper. 

In what follows, I consider different driving forces (SO, 
Zeeman, and hyperfine mediated), and the hyperfine spin 
relaxation mechanism that usually dominates in quan- 
tum dots, and derive mean field theory equations. It 
turns out that in all three cases the mean field theory re- 
duces to averaging over the longitudinal and transverse 
fluctuations of nuclear magnetization, in agreement with 
intuitive arguments. The results provide a justification 
for previous theoretical work performed in the framework 
of this approach. Mathematically, the derivation of mean 
field theory is based on relating separate terms of the 
power series for the electron spin flip probability W{t) to 
Eulerian F-functions of integer or half-integer arguments 
and employing their integral representations. The theory 
allows to derive, in the framework of a unified technique, 
time dependence W{t) of EDSR that is oscillatory for 
the SO and Zeeman mechanisms, in good agreement with 
the theory of Koppens et al.,^ and shows a nearly mono- 
tonic increase for the hyperfine mechanism in accordance 
with the conclusions by Laird et al^ In the latter case, 
W{t) approaches its saturation value Woo very fast, like 
a Gaussian exponent, according to the mean field theory, 
while corrections to this theory show a slow power-law 
decay. 

In this paper, the theory is developed as applied to a 
single quantum dot; second dot can only serve for prob- 
ing spin dynamics. The model includes neither electric 
current across the dot nor the coupling to thermal bath 
that can result in electrical pumping of nuclear spin po- 
larization and the back action of this magnetization on 
the electron spin.'^'^''^^ In electrically pumped double dots 
this back action is known to result in strong Overhauser 
fields and instabilities that were observed in the Pauli 
blockade regimej^s^'^'^^ For double dots, there is recent 
progress in the controllable electrical generating nuclear 
spin polarization22i^ and suppression of electron spin 
dephasing caused by random fluctuations ,-3§. With a dy- 
namically pumped nuclear spin polarization essentially 
exceeding the polarization from Poissonian fluctuations 
employed in Ref. ^ (about 50 G), one can expect dra- 
matic enhancement of the hyperfine-mediated EDSR. I 
expect that the techniques developed in this paper can 
be extended to those more general regimes. 

The paper is organized as follows. In Sec.|lTl the EDSR 
Hamiltonians are derived for three mechanisms of the 
coupling of electron spin to the driving electric field, and 
afterwards transformed to a form convenient for calcu- 
lating EDSR probability W{t). A general mean field 
expression for W{t) is derived in Sec. IIIII In Sec. llVi 
the equation for W{t) is simplified as applied to the SO 
and Zeeman (magnetic) mechanisms of EDSR, and the 
asymptotic behavior of W{t) is found. In Sec. |Vl the 
same program is performed for a more challenging prob- 
lem of hyperfine-mediated EDSR. Sec. IVII recapitulates 



the challenges and limitations of the mean field theory as 
applied to EDSR, and summarizes the basic approaches 
and results. Nuclear spin relaxation rate in a transient 
regime is calculated in Appendix A, where absence of its 
resonance enhancement is shown. Appendix B includes 
estimates of the corrections to mean field theory. 



II. THE HAMILTONIAN AND ITS 
TRANSFORMATION 



The Hamiltonian that is used in what follows is 
H = Ho + H,i{t) +Hz + HsoiKr) + Hhf{v) . (1) 
Here Hq is the zero-order Hamiltonian 



h k^ muj. 



0„2 



^0 = ^^ + ^r^ + V{y) , 
2m 2 



(2) 



describing an electron moving in (x, z) plane in a 
parabolic quantum dot, r — (x, z), with the in-plane con- 
finement frequency wq, being rigidly confined in y direc- 
tion by the potential V{y). Second term 

H^i (t) = eE(t) • r , E(i) = 2E cos cot , (3) 

describes the potential energy of the electron, with a 
charge (— e), in an in-plane driving electric field E(t) _L y. 
Third term 



Hz = -WbCB ■ s) , s = (t/2, 



(4) 



is the Zeeman energy for an electron with a negative g- 
factor, 5 < 0, like in GaAs and InAs, and s = a/2 and 
(T are the electron spin operator and vector of Pauli ma- 
trices, respectively, ^b = eh/2moc being the Bohr mag- 
neton. With a homogeneous field B in the confinement 
plane, B _L y, as in Refs.[3and[^, the diamagnetic contri- 
bution to the first term in Hq can be disregarded, hence, 
k will be identified as a canonical momentum. 

Generalized SO Hamiltonian i^so (k,r) consists of two 
terms, i?so(k,r) = Hsoi^) + Hso{r). Here i/so(k) 
is the usual momentum-dependent SO interaction, and 
only the linear in k terms will be considered in what 
follows. In the geometries of Refs. and [3, the field 
E(i) was applied along the face-diagonal direction, in our 
notations along [1,0,1]. In the related coordinate frame, 
the Rashba and Dresselhaus contributions to Hso can 
be written as Hji = afl(cr x k)y and Hd = au{<yxkz + 
o'zkx)- With z axis chosen along E(i), E = Ez, only the 
kz component of the momentum k matters in iJso(k). 
Hence, HsoO^) reduces to a single term 



HsoC^) = aa^kz, (5) 

'"^ Also, the homogeneous field B 



with a — an — o.r 

will be chosen as B = Bi because this is the simplest 

geometry in which Hsoi}^) of Eq. ^ mediates spin flips. 



The term _ff5o(r) describes SO interaction originating 
from the inhomogeneity of magnetic field, the mechanism 
employed in Ref. ^ With the polarization of the driving 
field E = Ez, only the z-dependence of Hso{^) matters, 
and to simplify calculations, Hsoi'''^) will be chosen linear 
in z. Keeping only the y component of the stray fieldj^ 



one arrives at 



Hso{r) = PcTy 



(6) 



where /3 is a slanting coeHicient. 

The standard expression for the hyperfine Hamiltonian 
is 



7I,./(r)-A^5(r-r,)(I,.s), A 



167r 
"3/ 



■W-BMn : 



where summation is performed over all lattice sites j, Ij 
are operators of nuclear momenta, / = 3/2 for GaAs, /i„ 
are magnetic moments of nuclei (difference in /i„ values 
for different nuclei is disregarded), and r/ is the enhance- 
ment factor. For GaAs, a rough estimate Auq w 10^'* 
eV can be used, with no = 4.5 x 10~^^ cm^^ for the 
concentration of nuclei,^ 



coupling and relates the EDSR intensity to it. If the spin 
resonance frequency ujg — \g\^BB/h is small compared 
with Wo, as is typical of GaAs, then lo/luq ~ '-^s/^o ^ 1- 
This is the special form in which the EDSR suppression 
due to the Kramers theorem— li^iiiiiS manifests itself as 
applied to parabolic quantum dots. 

Second, the canonical transformation of Eq. ^ 
changes the operator r as 



gik-R(t)j.g-ik-R(t) ^ J. 



R(i) 



(10) 



This time-dependent shift of r by R(t) transforms the 
r-dependent operators if 50 (r) and Hhf (r) as Hso (r) -^ 
Hso{r + R{,t)) and Hhf{r) -^ Hhf{r + R{t)). As applied 



(7) to 7fso(r) of Eq. ([6|), the time independent term pay 



produces EDSR only in conjunction with Hf,i(t), hence, 
this contribution is suppressed by the nonadiabaticity 
factor uj/ujQ and will be omitted. The dominant term 
in the transformed Hso (r) comes from the z-component 
of R(t) and equals 



H^{t) ^ payZ{t) . 



(11) 



It is not subject to the Kramers suppression because of 
the breaking time inversion symmetry. 



A. Transformation into the moving-dot frame 

In what follows, Hso{}^^ r) and Hfif{r) will be consid- 
ered as small compared with huo- However, before ap- 
plying perturbation theory, it is convenient to eliminate 
the zero mode inherent in the problem of a parabolic 
dot in a homogeneous field E(i). This mode manifests 
itself in keeping the shape of the electron cloud un- 
changed when it is displaced by an external homoge- 
neous and time-independent electric field4^ This can be 
conveniently achieved by performing a time-dependent 
canonical transformation of the Schrocdingcr equation 
ihdt^ = iJ* as 



2 . 
' 



^ir,y,t) ^ e-*-^WM'(r,y,t),R(i) = -eE{t)/mcu, 

(8) 

it describes changing to a coordinate frame moving with 
the dot. This choice of R(i) allows eliminating the driv- 
ing term Hei(t) and the zero mode. The transformation 
of Eq. ^ also produces a term —e'^E'^{t)/2mujQ that has 
no physical significance and can be eliminated by an addi- 
tional canonical transformation. However, two different 
results of the transformation to the moving-dot frame 
have important consequences. 

First, instead of the term Hei{t) a new term Hei{t) 
appears in the Hamiltonian 



Hei{t) = -2-^-^(k-E)sinwt. 



(9) 



This term, in conjunction with _ff5o(k), drives the SO 
mediated EDSR. A factor lo/ujo ^ 1 in Eq. ^ empha- 
sizes a critical role of nonadiabaticity for this type of SO 



B. Projecting onto the ground state 

[sec: project] At this moment, it is convenient to project 
the moving frame Hamiltonian onto the oscillator-type 
ground state of the zero-order Hamiltonian Hq. Two 
terms in the total Hamiltonian H, HsoO^) and Hei{t), 
are nondiagonal in oscillator quantum numbers. Pro- 
jecting them onto the oscillator ground state requires 
performing a standard (Luttinger-Kohn^ or Schrieffer- 
Wolff'*"' type) canonical transformation^ 



{Ho + Hi] 



\Ho+Hi)e- 



(12) 



with Hi = iJso(k) + Hei{t). Choosing T from the 
condition of canceling the linear in Hi term. Hi + 
[T,Hq] = 0, the Hamiltonian reduces to (_ffo + -ffi) — ^ 
Hq -— [T, [T, iJo]]/2 in the quadratic in T approxima- 
tion. Solving the equation for T results in (0|T|1) = 
— (O|i/i|l)/?ia;o, where |0) and |1) are standing for the os- 
cillator ground and first excited state, respectively. Then, 
using Eqs. ([5]) and ^ and keeping in the second order 
correction only the term oscillating at the frequency a;, 
one arrives at a Hamiltonian Hq + Hsoit) with 



Hsoit) = 2 



aeE Lo 



(Tx sin tot 



^^Erl 



iso Wo 



(Tj^sinwi. (13) 



This operator describes the joint effect of Hgo^) and 
Hei{t). In the projected Hamiltonian Hq = hujQ/2; this 
constant having no effect on spin dynamics will be omit- 
ted. When deriving Eq. P^ . the expression KOj/c^ll)!^ = 
mijjQ/2Ti for the matrix element of the momentum was 



used. In this equation tq = ^Jfi/mujQ is the electron 
ground-state radius, and (.30 = ^ /ma is a characteristic 
SO length. Eq. p^ holds when ea,eEro{(jj/uio) <C ?iwoi 
with £„ = ma^ /'n for the characteristic SO energy. 

Averaging _ff/i/(r-|-R(i)) over the ground state ^o(r, y) 
of the dot, expanding the average in R(i) and keeping two 
leading terms of the expansion, results in two hyperfine 
contributions to the Hamiltonian 

3 
Husit) - AY^mt) ■'^r,)^l{^„y,){l, ■ s): (14) 



here and below ^PoCr, y) is chosen real. The criterion of 
this expansion, Z{t) ^ rg, is equivalent to e-Erg ^ Uujq. 

The term H'^^ is time independent and describes the 
hyperfine corrections to the electron Zeeman splitting in 
the field B — Bz and the electron-mediated coupling 
between nuclear spins. It will be discussed in Sec. Ill Cl in 
more detail. 

The term Hhf{t) oscillates at the frequency oj and re- 
sults in hyperfine-mediated EDSR. Below, only step-up 
and step-down spin operators responsible for these tran- 
sitions will be retained in Hhf{t), hence, it reduces to the 
form 



HUt) 



Y, Z{t) ■ d.^ [^l{v,,y,)]{I+a^ + I-a+) , 

(15) 



± 



where a± = ax i: i<Jy, and If — H zt il^ 

Finally, the Hamiltonian in the quantum-dot reference 
frame reads 



H - [Hz + Hj-^f] 



+ [Hso{t)+HUt)+Hftf{t)]. (16) 



Three terms in second bracket represent three EDSR 
mechanisms discussed above. In deriving this Hamilto- 
nian, only linear in R(t) terms were retained. Terms of 
the higher order in R(i), as well as higher harmonics in 
Eq. (J13p and some different terms, require a theory in- 
cluding parametric excitations^^ 



C. Spin-spin coupling operator 

In this paper, electron spin coupling to the nuclear spin 
bath and the spectral diffusion caused by this coupling 
are considered as the only source of electron spin deco- 
herence because of the general scope of the paper focused 
on the hyperfine coupling and especially on the interplay 
between its contributions to different terms of the pertur- 
bation theory. This decoherence mechanism was recently 
discussed in the context of free spin induction decay and 
Hahn echo by Witzel and Das Sarma^^ and Yao, Liu, and 
Shamj^ see also references in these papers to extensive 
literature on the subject. Electron spin decoherence orig- 
inates from the spin dynamics in the nuclear bath caused 



by both the intrinsic mechanisms (like dipole-dipole in- 
teraction) and the extrinsic mechanism mediated by the 
hyperfine coupling of nuclear spins to the electron spin 
through the operator Hf^r of Eq. (|14p . In what follows, 
only the latter mechanism will be considered because it 
facilitates effective nuclear spin-spin coupling at the large 
spatial scale of vq. 



The term H^^ 



H, 



hf 



^H/^i includes two contributions. 



The longitudinal (secular) part ijj| r is proportional to Sz 
and describes random fluctuations of the electron Zeeman 
energy because of the Overhauser field. The transverse 
part H/[f, while depending on s±, cannot produce real 
electron spin fiip transitions because of the large electron 
Zeeman energy TilOs = I^Ims^ (nuclear Zeeman energy is 
small and will be omitted). However, in the second or- 
der of the perturbation theory in loJ^ , the operator H/^r 
results in nuclear spin nonconservation and spectral diffu- 
sion of the Overhauser field. Everywhere below, nuclear 
spins Ij will be considered as classical variables, i.e., their 
commutators will be disregarded. This simplifies calcu- 
lations and should not affect basic results because Ij and 
Ij' commute for j ^ j', and correlators involving more 
than two spin operators at the same lattice site are sta- 
tistically insignificant in the mean field theory developed 
below (Appendix B is the only exclusion). 

Applying the procedure of Eq. p^ to the operator 
Hz + H;^^, one finds T = Y.j aj(^>- ^ IJ (T+)/^/miJ:, 
with 



A 



y/2hL 



= *o(i-j>l/j), 



(17) 



and the transformed operator H/^r becomes an operator 
of the effective nuclear spin-spin coupling 



H, 



2r 



^n^m^n -^m^z • 



(18) 



It interchanges the projections of nuclear spins at differ- 
ent lattice sites while keeping electron spin unchanged. 
This expression for spin-spin coupling recovers the result 
by Yao et al^ Because the diagonal term, n = ?ti, is 
statistically insignificant for large dots, the n ^ m con- 
straint can be omitted. Then 



Hss = -^ g^ 



9-<^z , g± 



-E«.^ 



(19) 



In similar notations 

^1/ = 9z(Tz , g, = ^yhujs/2 Y «J^J ■ (20) 



We notice that g^ and g± have different dimensions. 

Substituting (fTOl) and (P0)l into ^TE\i , one arrives at the 
final form of the Hamiltonian in the dot frame 



H 



H, 



[Hz+Hlf 

[Hso{t)+Hra{t)+Htf{t)] 



(21) 



Three terms in the first bracket describe the electron 
Zeeman energy, its random shift due to the longitudinal 
component of the Overhauser field, and coupling between 
nuclear spins, respectively. Second bracket describes the 
driving force acting on the electron spin. It consists of the 
spin-orbit, inhomogeneous magnetic field, and hyperfine 
contributions. 

In conclusion of this section, one comment regarding 
the spin-spin Hamiltonian Hgs of Eq. (fT8)) should be 
made. It was derived by transforming the static Hamil- 
tonian Hz + H^p and the products anttm in -ffss are pro- 
portional to l/LOs^ with Ws playing a role of a large param- 
eter (spin gap) in the energy spectrum. However, in the 
EDSR regime the situation changes drastically because 
in the rotating frame Hamiltonian the large frequency los 
is reduced to [uJs — w), see Eq. ([M)) below. Therefore, 
a question arises whether Hss can experience a resonant 
enhancement at w — > Ws because of the "rotating frame 
singularity" . This problem is considered in Appendix A, 
where it is shown that Hss changes and becomes time 
dependent at the scale of the Rabi frequency but does 
not experience any resonant enhancement. Because solv- 
ing electron spin dynamics with a time dependent Hss (t) 
is an extremely challenging (or even impossible) task, in 
what follows the Hamiltonian Hss of Eq. (fT9|) is used 
as a model Hamiltonian. The magnitude of the coef- 
ficients Qj can be subject to renormalization, but their 
j-dependence will be chosen according to Eq. pT| . 



D. Rotating frame Hamiltonian 

All terms in second bracket of Eq. (PTjl depend on time 
harmonically, Hso{t) as sinwt while Hm{t) and H^Jt) 
as cosut. Applying the standard transformation to the 
rotating frame'^'^'^^ 

*(r, y, t) = exp {i(j,ojt/2)^ RF{r, y, t) , (22) 

that transforms Pauli matrices as 

a± (t) = a± exp {^iut) , (j^ (t) = cFz , (23) 

results in the rotating frame Hamiltonian 

H^ (^^s + 9z-\9+9-\. + \{I-CJ+ + Ua.). (24) 



ii^-E^.^ 



h = -4^.d^A<i^^^yi)] 



When deriving (|24p . the rotating wave approximation 
was applied with all fast oscillating terms omitted. Here 
As — Tiiuj — LOs) is detuning, gz is the Overhauser shift, 
and 5-|_(7_/2 describes second order coupling to the spin 
bath. 

Using ITTI) . P^ . and p^ . one finds driving terms /± 
for different mechanisms of EDSR 



/|^ = ±i(e^r^/^so)(^/^o), 



f^ = TiPi^E/mu^l) 



(25) 



(26) 



2tocJq 



(27) 



It is seen from Eqs. ([25]) - (P7)) that for two first mech- 
anisms the expressions for f± are very similar, hence, a 
unified theory of EDSR will be developed for them, see 
Sec. IIVI As distinct from them, presence of spin angu- 
lar momenta /■ in the coefficients f^ for the hyperfine 
mediated EDSR changes the situation drastically. This 
mechanism requires a special consideration, and a theory 
for this type of EDSR is described in Sec. |Vl 

Eq. ^^ for the Hamiltonian H, equations dUl) - ^^ 
for {gz^g±)^ and ([25| - ([27| for f± form the basic system 
of equations. In what follows, they will be solved and 
spin flip probabilities will be averaged over nuclear spin 
configurations for calculating time dependence of EDSR. 



III. SPIN FLIP PROBABILITY 

With the external magnetic field B strong enough, the 
initial electron spin state can be chosen as | t)i and 
electron wave function in the rotating frame evolves as 
^flp(i) = exp (-iff i/n) 1*0 T)- Then, applying ^ re- 
sults in the spin flip matrix element 

(i *o| exp {iazM/2) exp {-iHt/ri)\^Q t> 
and spin flip probability 

W{t) = |(i *o| exp {~iHt/h)\mo ])? , (28) 

because the factor exp {iazUjt/2) cancels out. The square 
of the Hamiltonian H of Eq. 

'As 1 



H^ = W 



9z 



9+9-] +f+f- (29) 



does not depend on Pauli matrices. Hence, the exponen- 
tial factor in (P5)l can be simplifled as 

exp{-iHt/h) = cos{m/h) - i{H/n) sm{m/h) . (30) 

Finally, the transition probability equals 



W{t) 



*o 



H^ 



■ sin^ (Ht/h) 



*o 



(31) 



This is the celebrated Rabi formula with /+/- for the 
driving term and (A^ -I- 2gz) for detuning;^- The spe- 
cial feature of Eq. PT|) is presence of the dephasing term 
—g+g- /2 originating from the transverse part of the ran- 
dom nuclear magnetization. 

Spin oscillations described by W{t) are controlled by 
the competition between the quantities Ag, gz, g+g-, and 
/+/_, all of them being small compared with the Zee- 
man energy Hz that does not appear in (|3ip explicitly. 
Corrections to the initial state *o| T) due to the nuclear 
spin fluctuations should merely renormalize the basic pa- 
rameters, in inverse powers of Hz, similarly to the Bloch- 
Siegert corrections to the rotating wave approximation,— 
without changing the main pattern of oscillations. Hence, 
they will be disregarded in what follows. 



A. Averaging over nuclear angular momenta 



The probability W{t) depends on the nuclear spin po- 
larization through gz and g± of Eqs. (fT9|) - ([20ll and f^-' 
of (|27p . Because experimental data are typically taken by 
averaging over dozens of thousand pulses covering time 
spans exceeding the nuclear spin diffusion timejii^ii^ av- 
erage values of W{t) of Eq. (pij) over all nuclear spin 
configurations are of the principal interest. Calculating 
these average values is highly facilitated by the fact that 
the series for the ratio 



sm{nt/h) 

n 



k=0 



2k 



■. (nt/h) 

{2k + l)\ 



(32) 



includes only even powers of Ti. Therefore, each term of 
the power series for the averaged probability W{t) 



wit)^ E 



f\k+k' /W^'\2(fe+fe' + l) 



X (/+/- 



^^^^ i2k + l)\{2k' + l)\ 

( — + .9^ - 2^+9- 1 + /+/- 


k+k' , 

) (33) 

/ O.nuc 



is a polynomial in the nuclear angular momenta Ij . The 
subscript {0, nuc} indicates that both quantum averag- 
ing over ^0 and statistical averaging over nuclear mo- 
menta should be performed. Due to the weak interac- 
tion between nuclear momenta, only single-site correla- 
tors can be retained. Next simplification originates from 
the fact that for large quanum dots, containing about 
10^—10^ nuclear spins, pair correlators statistically dom- 
inate; the contribution from higher correlators is esti- 
mated in Appendix B. Finally, of all pair correlators only 
(/|/|) = /(/ + l)/3 and (/+/7) = 2/(7 + l)/3 do not 
vanish. This allows to separate the averaging over gz 
from the averaging over g± and f± . 



B. Averaging over longitudinal magnetization 

From Eqs. (fT7|) and ([20| follows an expression for gl 
averaged over nuclear spin configurations 

(gl) = j^A^noHl + l)JJdrdy ^^(r, y) . (34) 

With oscillator ground-state functions "00(2^) = 
exp (— a;^/2rQ)/-\/7ri/-^ro and V'o(-2) in ^ and z directions, 
and hard wall confinement ip{y) = \/2/d cos (ny/d) in y 
direction, '^Q{r,y) = ipo{x)'4!(y)ijjQ{z) and 



{gl) = AV2 , A - .^A2„o/(/+l)/87rr2c 



(35) 



By combinatorial arguments based on the multiplicity 
of possible pairings, expressing {g1^^) in terms of pair 
correlators results in 



(<?: 



2M\ 



(2M-1)!!(5. 



2\M 



(36) 



with (2M - 1)!! = 1 for M = 0. Then, from the relation 
between (2M — 1)!! and the Eulerian Gamma function 

(2M-l)!! = 2^^r(A/ + l/2)/V^, (37) 

and from the integral representation of V{M + 1/2) 



r(M + ^ ) = / dii^M g-t (3g^ 



valid for integer values of Af , a Gaussian distribution for 
Q, follows 



2M\ 



(3 



Pa(w) 



A\AF 



dw w ' pa{w) 



exp {-w^ / A'^) . 



(39) 



with a standard deviation (g^) — A^/2. 

Therefore, averaging ([55)1 over the longitudinal magne- 
tization gz can be performed as a Gaussian integration 
with the distribution function pa{w) of Eq. pQ]) . 



C. Averaging over transverse magnetization 



Averaging Eq. ((33)) over /■ is more involved because 
these operators appear both in g± and f±. Applying 
(fT7)l and (fTO|) results in an expression similar to 
and finally results in 



(9+9-) = 



A^npljl+l) 
ATrhciJsrQd 



(40) 



Average values of {g+g-)^ , in the pair correlator ap- 
proximation, calculated from combinatorial arguments, 
are 

{{g+g^r)=M\{g+g^r^M\E^^'. (41) 

From the integral representation 



M! = T{M + 1) = / dtt 
Jo 



*^e-' , 



(42) 



by changing the variable t — iP' j c? ^ one arrives at a 
Gaussian distribution 



j2„ „,2M, 



{{9+9^r) = d'nu'^^'p^iu) 



Ph(m) 



expi-u'/E'). (43) 



Not surprisingly, the distribution is two-dimensional. 
This reflects the two-dimensional nature of the transverse 
spin polarization / . Constants A and S are connected 
by a simple relation 



A = -Jhujsl^ E. 



(44) 



Condition of spin resonance, A ^ hujs, is tantamount to 

SV2A < 1 . (45) 

This condition of weak spin-spin coupling will be used 
below. 

Next step is averaging ((55|) over f^. It follows from 
(|27l) that 

ift'f'i') ^ h.ns.i>{^)' 

X J J dr dy [d,[^l{r,y)]\ , (46) 
and after performing all integrations 



(/f/-^> 



no/(/ + l) fAeE\ 
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Trd \ huiQ J 



n' 



(47) 



The same combinatorial arguments that led to ((4T|) - (|43|) 
result in a two-dimensional Gaussian distribution 



1 



Pniv) 



2 /n2\ 



nQ' 



■ exp {-V /fl 



(48) 



cross-products g+f_ and g-fjj. According to (fT7|) . 



The last sort of pair products appearing in (j33[) are 
OSS-products 
dm), and (EZI 

(5+/!;/) (X / I dvdy 5.[*4(i., 2^)] ^ , (49) 



hence, cross-products vanish after the integration. Van- 
ishing the cross-products allows separating the averag- 
ing over the products g+g- from averating over fjj f_ . 
Potentially, cross-products describe the effect of driving 
force on spin dynamics of the nuclear bath, but in the 
approximation accepted their average values vanish. 



D. Mean field equation 



It was shown in Sees. IIIIBI and IIII CI that averaging 
over nuclear spins reduces to three integrations over the 
distributions pa(w), p~.{u), and pn{v). For these inte- 
grations, the variables in ((33| should be substituted as 

9z ^ w, g+g- -^ V?, f+f- -^ v^ . 

Afterwards, both summations over k and k' in (|33p can be 
performed, and one recovers the Rabi formula. However, 
now its arguments are auxiliary variables u, v, and w, 
over which integration should be performed. Finally 

W{t) = dw pa{w) / d^v pn{v) / d^u ps{u) 



,(50) 



smM i/f %+w-iu2) +v^ j- 



X V 



^+W-iu2 



with pa(w^), Ps{u), and pn{v) defined by ((39)) . (|43)) . and 
(j48]), respectively. Eq. ([50|) generalizes the equation ap- 
plied in Ref. l3 to hyperfine-mediated EDSR by including 
nuclear spin dynamics described by the term u'^ /2 and 
the corresponding integration. 

The above derivation proves that as applied to EDSR 
the pair-correlator approximation is tantamount to a 
mean field theory described by Eq. ((50|) . It was de- 
rived for the hyperfine-mediated EDSR. For the spin- 
orbit or magnetically mediated EDSR the integration 
over V should be omitted and w^ should be substituted 
by ff'-'f^'-' or /™/™, respectively. Corrections to the 
mean field theory are estimated in Appendix B. 



IV. SPIN-ORBIT AND MAGNETICALLY 
MEDIATED EDSR 



In this case (|50|) simplifies and includes integration two 
variables, u and w, with u^ substituted as f^ ^ /^ = 
/+/_ with fl^fl" or /;"/" of Eqs. (USD and ^ for 
the EDSR mediated by the spin-orbit or inhomogeneous 
magnetic field mechanisms, respectively. The results are 
also applicable to the usual spin resonance driven by an 
ac magnetic field. It is convenient, by using the relation 
sin^ (j) = {1~ cos 20)/2, to split W[t) onto its asymptotic 
value Woo at t — > cxD and the time dependent part Wi (t) 



W{t) = Woo-Wiit), (51) 

and introduce dimensionles variables 

(: = w/f, i^uy2f,T = ft/h, (52) 

and parameters 

6^Aj2f, R^f/A, R,^2f/E'. (53) 

Then 

RR, /■- f°° e^p{~R\^-R,0 

and 



Wi{t) = 



dCe 



-fl^C 



d^e 



-Rii 



2V^ J- 

((5 + C-0' + l 



(55) 



Woo can be calculated in the limit of weak driv- 
ing force, i? <C 1, due to the criterion of Eq. ([45]): 
this regime is typical of most experiments. Indeed, 
the integrand includes a narrow Lorentzian (5-function 

(5(C + (5-C) ~ 7r"V[(C + '5-0^ + l] that allows perform- 
ing integration over (^ for _R ^ 1. Afterwards, integra- 
tion over ^ can be easily performed because the Gaus- 
sian exponent is nearly a constant due to the fact that 
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2R^5/Ri = (A,/A)(S72A) < 1 for any A,/A < 1 due 
to (gS]). Finally, 



W^ 



f 



2y^A 



exp 



4A2 



(56) 



When (|45p is satisfied, Woo does not depend on S. Gen- 
eral expression for W^o can be found from (p5)l - (p5|) by 
omitting factors cos (2sr) in integrands. 



which allows expressing Wi(r) in terms of a Fourier 
transformation. Afterwards, integration over C in (|55p 
can be cut into parts in such a way that, after simple 
transformations like changing signs of variables, the inte- 
grations are performed over positive values of both vari- 
ables. Finally, one arrives at 



A. Simplifying the expression for Wi{t) 



rW, 



.(2) 



Wiir) = W^^'ir) + W^r(T) + W^^'ir) (58) 



(3), 



To simplify integral (j55p . it is convenient to introduce 
instead of ^ a new variable 



s^ = iS + C-0^ + l, s>0.^ 



(57) with 



W'^ Jo JVc^+i 



ds 



exp(-i?iVs2- 1) 
sVi2— y 



cos (2sr) 



(59) 



w['\t) 



RR\ 



dCe^v[-R^{C. - 5f - RiC\ 



Vc" 



'^ , exp(_Ri\/s^ 
ds ■ 



s\/s^ 



cos (2sr) , 



(60) 



2\/t^ Jo J\ 



exp(-i?iVs2- 1) 
sx/i2~T 



cos (2sr) 



r 



(61) 



After introducing two auxiliary functions 

G±(C) ^ I du exp [-i^^M ± (5)2 ± i?iu] 



(62) 



G'±(oo) = 0, equations ([55]) - ([SO]) can be performed by 
parts. This transformation reduces repeated integrals to 
onefold integrals because ^ enters into internal integrals 
only through their limits. Finally 



wi'Hr) 



RRi 



ds 



G+(V^^^T)-G+(0) 



20F Jo s\/i2~~T 

exp (-i?i v/s2 - 1) cos (2st) , (63) 



Wi'\r) 



20Fio sV^^^T 

(64) 



W/(3)(^) = _:^ /■°" G40)^^_«^^^^ 



e "" " 'cos(2st). 

(65) 

Because G±(C) can be expressed in terms of the Erfc- 
function, integrals (p^ - (p5|) can be easily calculated 



for any given set of parameter values. However, a wider 
insight comes from calculating their asymptotic behavior 
at large r. 



B. Large r behavior of Wi(t) 

Large-r behavior of Wi{t) is controlled by the sin- 
gularities of the integrands of the Fourier integrals (15^ 
- (j65p . All integrands decay exponentially at s — > cx), 
and square-root behavior at s = 1 is their only singular- 
ity. Because the numerator of the integrand of W^ (t) 
vanishes at s = 1, the singularity is softened. There- 
fore, the large r behavior is controlled by Wl (r) and 

Wi (r) due to Vs^ — 1 in the denominators. For calcu- 
lating the leading terms of the expansions, all other fac- 
tors in the integrands can be takes at their s = 1 values, 
~T) -^ G_(0), exp(±i?iVs2~~T) ^ 1, and 
V2\/s — 1. The remaining integral reduces 



as G^{\/s 



sy/s 



1 



to Fresnel integrals and can be easily calculated 



ds ,^ , / TT /„ vr , ,„„, 

cos (2st) = W — cos 2t + - (66) 

1 y/7^^ ' ' V 2t V 4 ' 



The coefBcient G- (0) equals 



2R 

i?r'exp(-A2/4A2); 



(67) 



second part of ([57]) follows from (|^5)) . Remarkably, while 
exact equation is asymmetric in 5, the approximate ex- 
pression is symmetric in the detuning Ag. Finally, for 

t — !■ CX3 



Wi{t) 



1 
2A 



■ exp 



4A^ 



2/ TT 

?i 4 



(68) 



Eq. (|68p was derived in a strict i — > oo limit. Be- 
cause of the existence in the problem of a large param- 
eter i?i ^ 1,R,R~^, one cannot exclude the existence 
of some intermediate asymptotic valid for t up to some 
upper bound related to this parameter. To check this 
possibility, one must first take the limit Ri —f oo in (|55p . 
Then 



Wi{t) = -^ / cosh{2R^6Vs^-l) 

V^ Ji sVs'^ - 1 

X cxp[-R^{s^ -1 + S^)]cos{2st). 



(69) 



Applying an approach similar to that leading to ([66|) , one 
arrives at ([68| . Therefore, it holds in the whole region 
1 ^ r < CX3. This type of universal oscillatory behavior 
with 1/vi decay was first found in the i?i — > oo limit, 
experimentally and theoretically, by Koppens et al..^^ 



HYPERFINE-MEDIATED EDSR 



To evaluate the three- fold integral of Eq. ([50]) . it is 
convenient to split it into two terms similarly to (|5ip and 



J 



W^i(r) = 



OO /'OO 



^ r dri T^'e-^' I dC e 
V^r Jo Jo 



rCHi / 3 _ ' 



define dimensionless variables 



C = w/n , r; = w/r2 , C = wVsrj , r = m/?i (70) 



and parameters 



R = f7/A , i?i = 2r2/s2 , 5 = A J2rj . (71) 



Then 



W^ 



RRi 



dC e-'"^ / dC e-«i« 

) JO 



and this integral can be evaluated in the R ^ 1 limit by 
performing integration over ^. Using (|45p one arrives at 



M^oo ~ (7rf]/4A) exp (-A2/4A2) 



(73) 



that is similar to (|56|) . More general expression for Woo 
can be found from (j8ip by omitting cos (2pr) in the in- 
tegrand. 

Evaluating Wi(r) is the major challenge, and rewrit- 
ing it as a sum of three integrals similar to (|58p (each one 
including integrations only over positive values of all ar- 
guments) is the first step. Transformations are similar to 
the applied when deriving (|5^ - (pTj) , but for the reasons 
that will become clear in what follows it is convenient to 
choose a different parametrization. In terms of a variable 
/?2 = (5 + ^ _ ^)2^ /3 > 0, chosen instead of ^, Eq. ([50]) 
reads 



cos{2^/]WTv^T) 



dC e 



-fl^(C-5)^-i?iC /"S^.fllg COs(2V/32+,y2^) 



RRi 



di] rj-^e '' / dC e" 
^0 



^fl=(C-<5)^--RiC 



d/3e' 



d/3 e- 



/32-f772 



_^^g cos(2v/^^T^r) 
/32 + 772 



(74) 



Introducing functions G± {(,) of Eq. (|62p and performing 
integrations over ^ by parts allows eliminating integra- 
tions over j3. The procedure is similar to deriving (|63p - 
([65|. Finally 



Wi{t) 



RRi 



drj T] e 



3.-n' 



with 



dC 



e + n 



F(C)cos(2v/C^T^t), (75) 



F(C) = G+(C)e'«^^ - [G-iOe''^^ + G^{0)e'"-<]. (76) 



10 



Functions G± {() are defined as complementary to G± (C) 



G±{C)= du e^p[-R^{u±Sf±Riu] (77) 



and obey a relation 

g±(-C) = -Gt(C)- 



(78) 



Eq. ([75]) is the final expression for the time dependent 
part of the transition probability. 



A. Large r asymptotic of Wi{t) 

Eq. (|751) demonstrates the convenience of the 
parametrization employed in ((74|) . Indeed, square root 
in the argument of cos (2-^/^^ + tj'^t) can be conveniently 
eliminated by transforming to polar coordinates as ^ = 
pcos0, 77 = psincj). Then ([75]) reads 



Wi{t) 



RRi 



7r/2 



X I dp p e 
'0 



■'^F(pcos0) cos(2pT). (79) 



The internal integral in (|79p is a Fourier transforma- 
tion, and its asymptotic behavior at large r is controlled 
by the singularities of the integrand. It is analytical for 
< p < oo. When both R,Ri ^ 0, F{() decays ex- 
ponentially with C, hence, F{p cos (j)) shows exponential 
decay for (f) 7^ 7r/2. The factor exp (— p^ sin^ (p) decays 
exponentially for 7^ 0. Therefore, the whole integrand 
decays exponentially for arbitrary value of (j). To find 
out the behavior of F{Q near ^ = 0, it is convenient to 
transform (|76p to an equivalent but a more symmetric 
form 



no 



[G+(C)e-^^^-G_(C)e^^'^] 
2G_(0)cosh(i?iC). 



(80) 



It is seen from the comparison of ([78l) and ([80|) that F{C) 
is an even function of C, ^(— C) = ^(0) and because it 
is analytical at C = 0, all its odd derivatives vanish at 
C = 0. This observation allows extending the integration 
over p to the whole real axis, —00 < p < 00, and suggests 
that the integral vanishes faster than any finite power of 
T. Because this statement is valid for arbitrary 0, it is 
valid also for the integral over (j). Therefore, the decay of 
Wi{t) is of the exponential type, and this is apparently 
the most general statement that can be made. As distinct 
from Sec. IIVBI I see no way for finding the t —> 00 
behavior before performing expansion in inverse i?i. 

Expanding F{Q in R^ provides a more specific out- 
look onto the asymptotic behavior of Wi{t) of ((79)) . Be- 
cause equations are rather bulky, results will be first 
provided for the resonance regime, S — 0. The lead- 
ing term of the expansion of F{C) is F{p cos 4>) ~ 



(2/_Ri) exp(— i?^p^ cos^ 0), and by substituting it into 
([7^ one arrives at 



Wiir) 



R /"^/^ f°° 

= —j= / d4> sin (j) dp p 

V'^ Jo J~oo 

X exp[-p'^{sm'^ (p + R^ cos^ 0)] cos (2pT).(81) 



The Gaussian integral over p can be performed exactly, 
but the result is rater cumbersome. It simplifies es- 
sentially when terms of the order of 1 are disregarded 
compared with the term t^ 3> 1. Afterwards, the in- 
tegration over <f> can be simplified by choosing a proper 
parametrization. Remarkably, for R = 1 the integrations 
over p and (j) separate, and Wi{t) shows a somewhat spe- 
cial behavior as will be seen from some examples in what 
follows. 

For i? < 1, by choosing 



1/[1- (l-i?2)cos2 



as a new variable, one finds 

Rt^ rVfl^ 



I^i(r)«- 



2(1 - i?2)3/2 y^ 



du 



l-R^u 



^/^I~T 



(82) 



(83) 



Because exp(— ur^), in the asymptotic region of large 
T, decreases very fast with u, the main contribution to 
the integral comes from the close vicinity of the lower 
integration limit, u — 1. However, the integrand in (j83[) 
either diverges in this point, for R ^ 1, or vanishes, for 
R—\. As a result, the asymptotic behavior of the inte- 
gral is strongly influenced by the competition between r 
and \R - 1|. For t'^{1 - i?) > 1 



Wi{t) « j^= re" 



(84) 



In the experiments by Laird et al.^ the coefficient R was 
small, R « 10~2. Nevertheless, it is instructive to con- 
sider a wider region of R values. Eq. ([84]) is valid for 
t2(1 — i?) ^ 1, but in the close vicinity of i? = 1, for 
1 ^ T^ ^ |1 — i?|^^, an intermediate asymptotic holds 



Wi{t) 



r^e" 



(85) 



For i? > 1, the following equation holds instead of 



Wi{t) 



Rt' 



2(i?2 - 1)3/2 J^^^^ ^ 



du 



R^u-1 



e^"^ , (86) 



and integration near the lower limit, for t^{R — 1) ^ 1, 
results in 



Wiir) 



i?2 



2(i?2 



1 



(87) 



Despite the difference in the specific form of the asymp- 
totic expressions of ([55]) - ([57)) . they have three major 



11 



features in common. First, in all cases the decay is of ex- 
ponential type, close to Gaussian, in agreement with our 
general conclusion. Second, as distinct from the Rabi os- 
cillations typical of spin-orbit and magnetically mediated 
EDSR (see Sec. lIVp.^'^ smooth time dependence is the 
distinctive feature of hyperfine-mediated EDSR.— Rabi 
oscillations are washed out by averaging over the fluctu- 
ations of the transverse nuclear magnetization. Third, in 
all cases Wi(t) is negative. This suggests that the prob- 
ability W{t) of ([5T|) approaches its limit value, at i — ^ cx), 
from above, in agreement with a weak overshoot found 
numerically by Laird et al.^ 

Corrections to ([8T|) from the higher order terms of 
the expansions of G±(C) and G±(C) in \/Ri contain the 
same Gaussian exponent exp {—R^p^ cos^ 0). Their pre- 
exponential factors are small in [R/ Ri)^ = (S^/2A)^ and 
{R^C,/Ri)^ what is actually the same because RC, ^ 1 due 
to the Gaussian exponent in (j8ip . These corrections can 
also be estimated in terms of the time r, after perform- 
ing both integrations in ([5T|) . They are of the order of 
[R^T/Rif = [(17/A)(SV2A)]2t2 and can become im- 
portant only for very large values of r because the esti- 
mate includes, in addition to the small parameter S^/2A, 
Eq. ([45l) . also a factor 0/A that is typically small. 

Eq. ((8T|) was derived for exact resonance, (5 = 0. For 
(5 7^ 0, the leading term of the expansion of F{C,) in 
l/Ri < 1 acquires a factor exp {-5'^R^) cos {25R^C). Af- 
ter performing the Gaussian integration over p in the 
modified equation (|8ip . the leading term of the expan- 
sion in T^ ::^ 1 acquires a factor cos {25t/ cos 0) . After 
changing the variable from </> to u according to (|82p . a 



factor cos [25yju[l — R'^)/{u — 1) r] should be incorpo- 
rated into the integrand of Eq. ([55)1 . Finally, one arrives 
at a generalized equation ([M)) 



Wi{t) 



i? V^ exp [-(52(1 _ij2)] 



re 



(88) 



2VT^i? 

Remarkably, for i? ^ 1 the suppression of EDSR due 
to the factor e^^ is much stronger than in (p5)) where 
the suppression factor e~^ ^ . However, for R — 1 the 
suppression disappears at all. This is another manifesta- 
tion of the special EDSR regime near the R = 1 point 
where the driving force coincides with the characteristic 
spectrum width. 

One more difference between the off-resonance regime, 
6^0, and the on-resonance regime, (5 = 0, is in the 
structure of the 1/Ri expansion series. The expansion 
coefficient R^t/Ri is the same in both cases, but for 
6^0 the series includes all powers of R^t/Ri, odd and 
even, while in the resonance regime odd terms vanish. 

In all cases the characteristic time of EDSR, when 
W{t) reaches its flat maximum, is r ^ 1 or i ^ h/H.. 

VI. CONCLUSIONS 

As distinct from ferromagnets, where macroscopic 
magnetization is stabilized by strong exchange interac- 



tion, direct interaction of nuclear spins in quantum dots 
is very weak. Therefore, the description of nuclear po- 
larization as a macroscopic field acting on electron spin 
is based on a large number of nuclear spins inside the 
dot, iV 3> 1, and their slow dynamics compared with 
electron spin dynamics. Meanwhile, the hyperfine cou- 
pling of electron spin to nuclear spins is usually much 
stronger than the direct coupling between nuclear spins, 
and under the conditions of EDSR the Zeeman term van- 
ishes in the rotating frame. Therefore, applying macro- 
scopic approach to EDSR needs justification. Moreover, 
only consistent microscopic approach to EDSR can en- 
able calculating corrections to macroscopic theory and 
provide solid basis for including nuclear dynamics and 
investigating instabilities in electron-nuclear system. 

From the standpoint of the Central Limit Theorem, 
one anticipates that EDSR probability can be found by 
averaging the Rabi formula 



Wnit) = f 



sin^ W{e/2y + p t/h] 
(6/2)2 + /2 



(89) 



over the Gaussian distribution of random nuclear mag- 
netization. However, having in mind that both detuning 
e and driving force / depend on nuclear magnetization, 
performing the averaging of Wfl(t) is a challenging prob- 
lem per se. Deriving this recipe from microscopic theory 
and developing techniques for going beyond the macro- 
scopic description is even a much more challenging task. 

The approach developed in this paper is based on the 
following main steps. First, because the expansion of 
sin (Tit/h)/Ti. in power series includes only even powers of 
Ti, Eq. (|32[) . the quantum- mechanical expression for spin- 
flip probability W{t) can be expressed in terms of poly- 
nomials in components of nuclear angular momenta Ij. 
Second, retaining only pair correlators (an equivalent of 
mean field theory), the contributions of the longitudinal 
and transverse magnetization can be decoupled, and sta- 
tistical weights of different contributions can be expressed 
in terms of F-functions of half-integer and integer argu- 
ments for longitudinal and transverse magnetizations, re- 
spectively. Third, employing integral representations for 
these F-functions, one arrives at Eq. (j50p describing one- 
dimensional averaging over longitudinal fluctuations and 
two dimensional averaging over transverse ones. Fourth, 
higher correlators can be accounted for in a regular way, 
see Appendix B. 

The probability W{t) can be split into its saturation 
value Woo and time dependent part Wi{t), Eq. ([?T|) . 
For spin-orbit and inhomogeneous Zeeman-coupling me- 
diated EDSR, Wi (t) can be reduced to onefold integrals 
([63)1 - (|65)) ; its large-i behavior is described by oscillations 
showing a i"^/^ decay in agreement with Ref. [S^. For 
hyperfine mediated EDSR, time dependence is close to 
monotonic, with a weak overshoot as found in Ref.|8|, and 
is described by a twofold integral f79\i . Large-i asymp- 
totic is Gaussian in the framework of mean field theory, 
but it develops a power-law tail when the first four-spin 
correlator correction is accounted for, see Appendix B. 
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The contribution of such a tail should become essential 
at moderate times i ex IniV, where N is the number of 
nuclei in the dot. 

Regarding the effect of nuclear spin-spin coupling on 
EDSR, it is of critical importance that this coupling 
shows no resonant enhancement under the EDSR con- 
ditions, see Appendix A. As a result, it has no effect on 
EDSR as long as the criterion (P5|) is fulfilled. Its explicit 
form is 



"2 

2A 



SNhuJohujs 



N — Trr^driQ . 



(90) 



For GaAs quantum dots with Ano w 0.1 meV, huja w 0.1 
meV, d « 5 nm, and B « 100 G, this ratio is small, 
SV2A- 10-2 < 1. 

It is the main restriction of the model employed in this 
paper that it includes only "passive" processes in the 
nuclear subsystem like spin relaxation, and does not in- 
clude mechanisms that can result in pumping nuclear po- 
larization by electron dynamics. Including such "active" 
processes and related instabilities in the coupled electron- 
nuclear system into the formalism developed above is the 
major challenge. 
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APPENDIX A: DYNAMIC NUCLEAR SPIN-SPIN 
COUPLING 

Eq. (fTg]) was derived for a static Hamiltonian Hz + 
H^t. For exploring the role of driving force on spin-spin 
coupling, one can use a model Hamiltonian 



H 



-huj.a. 



1 

2 


Ihujs 

V 2 


3 


- + ^i^+) 


1 
2 


/(^+ 


e'-^^ + a.e- 


iujt\ 


(Al) 


de 


Qtical to Hz + H^j , 


with the 



second term rewritten in terms of the coefficients aj of 
Eq. P?|) . Last term describes the driving force in the ro- 
tating wave approximation, with / defined as / = /± of 
Eq. (j24p . The corresponding rotating frame Hamiltonian 



H = Hq + h{t) , Ho = -A,,cr^ + /ctx 



m = 



flUOs 



5]a,(/+a-e-*+/7a+e— *).(A2) 



For finding a dynamic analog of p^ . effect of h{t) on 
slow nuclear dynamics should be calculated in the second- 
order perturbation theory. To this end, by rotation in the 
spin space 

0-1 - 70-3 0-3 +70-1 /A ON 

(Tx = , „ , CTj^ = era , cr^ = , (Ad) 

\/l + 7 Vl+7 

with 7 = As/2/, the Hamiltonian H^ simplifies as 



Ho = fa,, /'-/VT+7^, 



(A4) 



and h{t) can be transformed into the interaction repre- 
sentation h{t) = exp {iHQt/h)h(t) exp {—iHot/h) by us- 
ing 



cri(i) = ai 



^3(0 



0^2 
0-3 



aa sin {2 ft) , 
a2 sin {2 ft) ; 



(A5) 



here and below, /' in frequency units. Then, calculating 
the second order correction through a i-ordered product 



til 
2! 



T{h{ti)h{t2)}dtidt2 



(A6) 



eliminating all fast-oscillating terms proportional to 
e^"^*, and exponentiating the result, one finds a hy- 
perfine correction to the phase of the electron-spin 
wave function. Representing the phase factor as 
exp[— ziJss(i)i/?i] allows finding effective spin-spin cou- 
pling Hss{t)- 

Exact expression for Hgg {t) is rather bulky, but in the 
large uj limit it essentially simplifies 



Hss{t) 



1 UJs/lU 



2\/l+7 
sin(4/'t) 



2 ^-^ 

vrin 



%0,nJ^rn^n 



Aft 



o-i 



sin^ {f t) 
ft 



-(72 + CTS 



(A7) 



Near the spin resonance ui ^ ujs, and the factor lUs/lo ^ 1 
can be omitted. Hss{t) oscillates with the Rabi frequency 
/' and saturates for ft -^ cxd. Expressing {<7i,a3,a^} 
through {ax,cfy,az\ by using (jASp . one can check that 
at short times, ft <C 1, Hss{t) coincides with Hss of 
Eq. din]). In the opposite hmit, for 7 > 1 and /'t > 1, 
Hss{t) ~ 5Hss/2. 

Therefore, Hss{t) shows no resonant enhancement at 
Lu — ^ LUs and has a magnitude comparable to H^s of 
Eq. (|19p and the same sign. In physical terms, this result 
suggests that the hypcrfinc dynamics near the resonance. 
As ~ 0, is strong enough to smear the rotating frame sin- 
gularity and reduce its effect to moderate changes in the 
nuclear spin relaxation rate that acquires oscillatory be- 
havior. Regarding the significance of this conclusion, see 
discussion at the end of Sec. Ill CI 
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APPENDIX B: EFFECT OF HIGHER 
CORRELATORS 

Eq. (j50p was derived including only pair correlators of 
nuclear spins (1^1^) = {IaIp)Smn, Sec. IIII Al In this 
Appendix, effect of higher correlators on EDSR inten- 
sity will be estimated as applied to the operators f^ of 



Eq. ([27|) . To simplify bulky calculations, terms gz and 
g+g^ in the Hamiltonian of Eq. ((24|) will be omitted. 

In the pair-correlator approximation, a correlator 
Fm = iif+f-)^^) including M pairs of operators f±, 
was evaluated as 



J^M = M!(/+/_ 



\M 



(E^") 



M 



(Bl) 



This estimate came from counting the number of possible 
^n-Jn pairings without any constraints imposed on m and 
n. However, such a count overestimates the contribution 
of pair-correlators that should be found from 



F^°^ = M! < /_/+ > 



M 



E 



"lZ lz lz 

mi 7712 " fnM 



(B2) 

Indeed, each sharing of indices like mi = TO2 results in 
higher correlators like ((/+/_ )'^). Therefore, one needs 
to evaluate the effect of this constraint on W{t). For 
this evaluation, it is enough to keep in (jB2[) only the 
principle sequence of corrections that includes terms with 
a number of pairs of shared indices, like mi = m2, but 
does not include triple shared indices, like mi = m2 = 
rris, and terms with a larger number of shared indices. 
Their contribution is of higher order in 1/-/V, where N = 
irr^driQ is the effective number of nuclei in the dot. 

Each term of the correlator J-j^j involving s pairs of 
shared indices TOo, makes a contribution 



-yj'M 



Ml 



s!2''(M-2s) 



X^h4 



VZ^ 



'/(E^^ 



2 \2s 



(B3) 



Here < s < Smax, with Smax = M/2 when M is even 
and Smax = {M — l)/2 when M is odd. The factor 
with factorials equals the multiplicity of independent se- 
lections of s pairs of shared indices mj. The origin of 
the sign factor {—Y becomes clear when one eliminates, 
in succession, first a single pair of shared indices, then 
the second, etc. Last factor can be calculated using the 
definition ^ of b^ with ^o(r,2/) of Sec lHIBl 



E^t/(E^" 



35/24iV. 



(B4) 



(0) 



Therefore, in the principle sequence approximation, Fj^j 
equals to 



M! / 35 \ 



' ^ ' hi s\{M-2s)\\A^n) 



(B5) 



Estimating a few first terms of the series shows that the 
expansion parameter is A'P /N . It corresponds to a typ- 
ical fluctuation M ^ ^/N in the nuclear reservoir. With 
M > 1 and s < Af , the ratio M\/{M - 2s)! « M^", and 
(IB5[) is an exponential series. It converges fast, and one 
can extend summation to infinity. Then 



M 



Tm exp (-35MV48Af) . 



(B6) 



Therefore, constraint (|B2p results in an exponential fac- 
tor of Eq. ((B6)) . 

For calculating the effect of the constraint (|B2p on 
W{t), Eqs. dHI]) and (|B6]) should be plugged into ^ 
(with gz,g± =0), and M\ should be eliminated by apply- 
ing averaging over the Gaussian distribution of Eq. 
Then one arrives at 



W^"\t) « /d^vpoH E 



fc,fc'=0 

h v\{k + k' -v)\ V 2 



(\k+k' /^ /^\2(fc+fc' + l) 

(2fc + l)!(2fc' + l)! 



2M / 35Af2 

V exp 

^' 48iV 



, (B7) 



M=k+k'~v+l 



where the superscript zero indicates that W^^' includes 
only pair correlators. With the last exponential fac- 
tor omitted, the sum over v \s a. binomial i;^[(As/2)'^ + 
^2jfc+/£ ^ g^j^j (|B7p reduces to a simplified version of ((50|) 
with only a single integration. 

In resonance, A^ — 0, the only characteristic time scale 
is fi/Vt. For T = t^/h ^ 1, sums over k converge fast, 
at k,k' ~ 1. Therefore, M ^ \ and the effect of the 
exponential factor exp(— 35M^/48iV) is small, only about 
1/iV. For large detuning, A^. 3> ri, a short time scale 
?i/As develops. However, it was shown in Sec. IV A] that 
EDSR frequency fluctuations A/?i of a scale A > As/2 
[that were omitted when deriving (jB7[) ] reestablish the 
fi/Q, scale, hence, the \/N estimate should remain valid 
for the hyperfine-mediated EDSR. 

The difference between F^^i and J^m reflected the re- 
duction in the contribution of two-spin correlators, for 
given A/, because of the constraint (|B2p . Next step is 
estimating the direct contribution of higher correlators, 
and this will be done for A^ = 0. The contribution of a 
single four-spin correlator ((/+/_)^) to W{t) equals 



iy(°)(t) « E 

fe,fe'=0 

7 

In 



(_)/c+fc'(^/^)2(fc+fe'H 



^. 



(2fc+l)!(2fc' 
M{M - 1)" 



1)! 



M- 



(B8) 



M=k+k' + l 



The coefficient is a product of the ratio 
((/_i_/_)^)/(/_i_/_)^ — 6/5 evaluated for classical spins 
and the ratio of 6-sums of Eq. (|B4p . The contribution 
of pair correlators was estimated in the principle se- 
quence approximation, sinfilarly to (|B3p . The multiplier 
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M{M — l)/4 reflects change in the combinatorial factor 
due to adding a single four-spin correlator and reducing 
the number of two-spin correlators by two (total number 
of spin operators kept fixed). Exponential factor about 
exp (-36MV48A^) is omitted. 

The factorial A/! in J-m can be eliminated by intro- 
ducing a Gaussian variable v and performing Gaussian 
integration according to (|15)) . Afterwards, the factor 
M{M — 1) can be eliminated in the integrand by sub- 
stituting it by derivatives like Mv^ = vd{v^^)/dv. This 
transformation allows performing summations over k and 
k' and one arrives at 



W'^^\t) 



7 



d{v^) 



dv v'e-'""!^" 



sm^ivtjK) 



d(w2) 



(B9) 



The integral can be simplified by performing by parts. 
Finally, in terms of r = tQ./h, 

VF^i) (r) = — - / du u{2 ~ Au^ + u'')e-"' sin^ ut . 

(BIO) 
Its asymptotic behavior for t ^ 1 



the power-law decay originates due to finite slope of the 
integrand at u = 0. Therefore, at the time span of r ^--^ 1, 
first correction to mean-field behavior of W{t) from the 
four-spin correlator is of the order of 1/A^, i.e., of the 
same order as the correction (jB7|) from the reduction of 
the number of pair correlators. A new feature is the time 
dependence of VF'^-'(t). While the Gaussian asymptotic 
of Eqs. ([M)l - ([57|) is a distinctive property of the mean- 
field theory of hyperfine-mediated EDSR, (jBlip indicates 
a much slower decay of the corrections to this theory. The 
power-law asymptotic deserves a more detailed study, 
however, a cancellation of power-law terms and return 
to the exponential asymptotic behavior seems improb- 
able. And power-law tails suggest deviations from the 
mean-field behavior at relatively modest times of about 
T-lnA^. 



W^(i)(T)»7/(167Vr2) 



(Bll) 



It deserves mentioning that (jBTj) has implications be- 
yond the simplified model of the hyperfine-mediated 
EDSR for which it was derived. Spin-orbit and mag- 
netically mediated EDSR's show oscillatory behavior, 
Sec. IIVBI and because the convergence of the series 
describing oscillations slows down rather fast with the 
number of oscillations, this immediately translates in 
the increase in the corrections that are proportional to 
M^ K, [k^k'Y according to JBTj. Therefore, the magni- 
tude of the corrections to mean field theory is expected 
to increase essentially with the number of oscillations. 
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